Short time aging in binary glasses 



Giorgio Parisi 
Dipartimento di Fisica, Universita La Sapienza 
INFN Sezione di Roma I 
Piazzale Aldo Moro, Roma 00187 

February 1, 2008 

Abstract 

We present some simple computer simulations that indicate that at short time aging 
is realized in a simple model of binary glasses. It is interesting to note that modest 
computer simulations are enough to evidenziate this effect. We also find indications of 
a dynamically growing correlation length. 

Aging has been discovered long time ago and it is has been experimentally studied in great 
details in some materials . It has been later realized that it is a quite common phenomenon 
in physics. It has been the subject of wide theoretical investigations only recently [0, El ffl. In 
the nutshell aging predicts that the response of the system to a force that has been applied 
for a time t depends on t (also for very large times) when t is comparable to the waiting time 
t w , i.e. the time at which the systems has been carried to conditions at which the experiment 
has been done. 

Aging has being studied analytically in generalized spin glasses: in the simplest form 
it predicts that the correlation functions among a configuration at time t w and at time 
t w + t depend only on the ration t/t w in the limit of large times ||. This form of aging 
is approximately found to be correct in spin glasses both in experiments and in numerical 
simulations || [| [7j, (although some small modification may be needed). Slightly different 
forms of aging have been proposed, e.g the scaling variable could be t/tf^ with fi near but 
not equal to 1. Here we stick to the t/t w scaling and we will refer to it by the name of simple 
aging. 

The aim of this note it to start a systematic study of aging using numerical simulations 
in real glasses. We will show that the aging regime starts at relatively short times and 
some of its properties can be investigated with a modest amount of computer time. We 
limit ourselves to the analysis in the initial time region, leaving to further more systematic 
investigations the study of the behaviour at larger times and in bigger systems. 

The numerical experiment we present here is rather simple: we run a numerical simulation 
where the system starts from a fully random configuration (i.e. at infinite temperature). The 
system is then carried (at time zero) at temperature T. We make a photograph of the system 
at time t w and we compare the later evolution of the system with this reference configuration. 

The main quantity on which we concentrate our attention is the two times correlation 
function 
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Figure 1: The energy density for 258 particles as function of time. The fit is E — + At 



with E^ = 1.92, A = 12 and A =.7. 



where N is the total number of particles and Xi(t) denotes the position of the particle % at 
time t. In the limit of large times g(r,t,t) goes to the usual correlation function for liquids 
(apart from an extra delta function at the origin, which is absent in the usual definitions 
where the sum is restricted to i ^ k). 

For this function simple aging predicts that for both t and t w large, we have the scaling 
relation 

g(r,t w ,t) G(r,s), where s = — . (2) 

We have tried to test this relation for binary fluids. The model we consider is the 
following. We have taken a mixture of soft particle of different sizes. Half of the particles 
are of type A, half of type B and the interaction among the particle is given by 



(3) 



where the radius (a) depends on the type of particles. This model has been carefully studied 
studied in the past g, ||, [UJ. It is known that a choice of the radius such that ctb/cta = 1.2 
strongly inhibits crystallisation and the systems goes into a glassy phase when it is cooled. 
Using the same conventions of the previous investigators we consider particles of average 
diameter 1, more precisely we set 
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(4) 



Due to the simple scaling behaviour of the potential, the thermodynamic quantities 
depend only on the quantity T 4 /p, T and p being respectively the temperature and the 
density. For defmiteness we have taken p — 1. The model as been widely studied especially 



for this choice of the parameters. It is usual to introduce the quantity T = f3 4 . The glass 
transition is known to happen around T = 1.45 0. 

Our simulation are done using a Monte Carlo algorithm, which is more easy to deal with 
than molecular dynamics, if we change the temperature in an abrupt way. Each particle is 
shifted by a random amount at each step, and the size of the shift is fixed by the condition 
that the average acceptance rate of the proposal change is about .5. Particles are placed in 
a cubic box with periodic boundary conditions and at the end of each Monte Carlo sweep all 
the particles are shifted of the same vector in order to keep the center of mass fixed [11|. This 
last step in introduced in order to avoid drifting of the center of mass and it would be not 
necessary in molecular dynamics if we start from a configuration at zero total momentum. 

In our simulations we have considered a relatively small number of particles, N = 34, 
N = 66 and N = 258 (most of the data we show are for N = 258). We start by placing the 
particles at random and we quench the system by putting it at T = 1.8, i.e. at a temperature 
well below the glass transition. The energy as function of the Monte Carlo time (t) is shown 
in fig. 1. The data are the average over 25 different realization of the dynamics with different 
initial conditions. The energy seems to decay to an asymptotic value with some corrections 
which vanishes as a power of time. 
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Figure 2: The correlation function G(r, t w , t) as function of r for N = 258 for different values 
of t w (t w =128, 512,2048 ,8192, at t/t w equal to 3. The data have been averaged over 25 
different initial conditions 

We have measured the correlation function G for different choices of t w (i.e. t w =32, 128, 
512, 2048, 8192) at s = t/t w equal to 3. If we exclude the points at t w = 32 (which we have 
not plotted), the data for the correlation functions have a very similar shape independently 
of t w at fixed s, as expected from simple aging. In order to check in a more quantitative way 
the aging and to evidenziate possible violations of aging, we have introduced the quantity 
q(t w , t) defined as 



q(t w ,t) = J dxg(x,t w ,t)f(x) = 



i,k=l,N 



f(Xj(t + t w ) -x k (t w )) 
N 2 



(5) 



where we have chosen the function / in such a way that it is sensitive to the area of the 
central peak, i.e. 

f{x) = x^^> (6) 
with a = .3 The function / is very small when x » .3 and near to 1 for x < .3. The value 
of q will thus be a number very near to 1 for similar configurations (in which the particles 
have moved of less than a) and it will be much smaller value (less than .1) for unrelated 



configurations; using the same terminology as in spin glasses |P2| , |T5| , [Tj] q can be called the 
overlap of the two configurations. 
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Figure 3: The overlap q as function of s for t w =128, 512, 204£ 
258, averaged over 25 initial conditions. 
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In fig. 3 we plot the overlap as function of s at different value of the waiting time. We 
notice that for s near to zero there is noticeable dependence on the waiting time, which can 
be related to the fact that t is not large. There is a small drift of the data at large s. In 
order to see evidenziate the effect we plot in fig. (4) the function q at s = 4 as function of t w . 

It is not clear if the increase at larger values of t w is due to finite size correction or to 
small violations of the simple aging hypothesis. This point should be carefully investigated 
in the future and its clarification goes beyond the aim of this note. The data on smaller 
systems indicates that there are definitely strong volume effects at large times fTS] . This 
may seem surprising; in order to find the physical origine of this effect it is instructive to 
look to the time dependance of the fluctuations of q from sample to sample, i.e. to 

A(t w ,t) =< (q- <q>) 2 > 



(7) 



where the average is done over different realization of the initial conditions. This quantity 
seems to increase a a power of t w with and exponent \i near to 1/2. In fig. 5 we show the 
data (at s = 4) for 



M(t w ) = NA(t W} 4t v 
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Figure 4: The value of the overlap q at s = 4 as function of t w 
for iV = 66 and iV = 258. 

versus t]/ 2 . The quantity M(t w ) seem to increase a a power of time and shows a weak 
dependance on the size (as expected, usually fluctuations goes to zero as iV -1 / 2 ). 

In order to interprete this results it is convenient to recall the physical picture at the 
base of simple aging: the system evolves by a sequence of quasi equilibrium states and 
remains in a given state for a time proportional to time needed to arrive to it. In the most 
extreme picture we have a punctuated equilibrium of long period of stasis intermixed by fast, 
thermally activated, tunnelling events. Increasing the value of t w the barrier that we have 
to cross become higher ad higher and such collective movements involves a large number of 
particles. Roughly speaking we expect that the variance A is inversely proportional to the 
number (N/M(t w )) of regions which have been moved independently, i.e. 
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The previous result imply that the time variation of q is dominated by events which 
involve the rearrangements of regions of size which increase at least as W 3 . This type of 
behaviour (i.e. a dynamical correlation length increasing as a power of time) has been seen 
in quenched disordered systems like spin glasses [[U| |T7p . It is quite evident that there is a 
change in the behaviour when the number particles involved in a typical rearrangement be- 
comes of the same order of the sample size. Independently from any theoretical speculations 
it is clear than when the variance A becomes comparable to (1 — q) 2 the distribution of q 
cannot be anymore Gaussian (the overlap cannot become greater than 1!) and we enter in 
a new regime. Strong finite size effects are thus expected for sufficient large time, i.e. for 
times which increase as N 1 ^. It is also quite likely that for sufficient large times a finite 
system reaches a sufficient low state such that further jumps are inhibited or happen on a 
much lager time scale 



Summarizing the conclusion are this note are the following: 

• Aging effects can be observed also at relatively short times in binary glasses. 
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Figure 5: The value of M{t w ) as function of t l J 2 at s = 4. 



100 



Simple aging is observed for a variation of about two order of magnitude in the waiting 
time with deviations with are at most 10%. 

It is suggested that there is a dynamical correlation length, that indicate the size of the 
regions which are collective rearranged, which diverges as a power of time. (A direct 
study of the size of the rearranged regions in the equilibrium dynamics can be found 
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Further numerical simulations are needed to decide if the small violations of simple scaling 
are finite volume effects or they survive in larger samples. It would be also interesting to study 
the temperature dependence of the effect and compare the results with detailed theoretical 
predictions. 
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